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Abstract. A tomographic method is described that makes it possible to follow the propagation of shock waves 
across the photosphere of long-period variable stars. The method relies on the correlation of the observed spec- 
trum with numerical masks probing layers of different atmospheric depths. The formation depth of spectral lines 
is derived from synthetic spectra of non-variable red giant stars. When applied to Mira stars around maximum 
light, the tomographic masks reveal that the deepest photospheric layers are generally characterized by blueshifted 
absorption lines (attesting their upward motion), whereas the uppermost layers generally exhibit redshifted ab- 
sorption lines (indicating their infalling motion). Double absorption lines are found in intermediate layers, where 
the shock front is located. At later phases, the shock front is seen moving towards upper layers, until it leaves the 
photosphere. 
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1. Introduction 

Besides the brightness fluctuations, long-period variable 
stars (LPVs) also exhibit striking spectral changes, like 
the doubling of several absorption lines around maxi- 
mum light, first reported by Adams (1941). The line- 
doubling phenomenon has been source of conflicting theo- 
ries. Schwarzschild (1952) has suggested (originally in re- 
lation with W Vir Cepheids) that the doubling of the ab- 
sorption lines around maximum light is related to the pas- 
sage of a shock wave through the photosphere. Alternative 
models accounting for line doubling without resorting to 
differential atmospheric motions have been proposed by 
Karp (1975) and Gillet et al. (1985). 

In Paper I of this series (Alvarez et al. 2000), it was 
shown that the temporal evolution of the red and blue 
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peaks of the double absorption lines of the Mira vari- 
able RT Cyg follows the "Schwarzschild scenario" . This 
result thus indicates that the line-doubling phenomenon 
is caused by a shock wave propagating in the photosphere. 
In Paper I, it was shown that this propagation may even 
be followed by resorting to a tomographic technique. The 
method relies on the correlation of the observed spectrum 
with numerical masks probing layers of different atmo- 
spheric depths. This technique applied to RT Cyg revealed 
that the evolution of the line doubling with atmospheric 
depth is also consistent with the Schwarzschild scenario. 

This paper presents the implementation of this tomo- 
graphic technique. Its application to a large sample of LPV 
stars is presented in Paper III. The possible relationship 
between absorption line-doubling and emission lines will 
be investigated in a forthcoming paper in this series. 

2. Scanning different parts of the atmosphere 
with the default ELODIE correlation masks 

A monitoring of a large sample of LPV stars was per- 
formed with the fibre-fed echelle spectrograph ELODIE 
(Baranne et al. 1996) at the Observatoire de Haute- 
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Vr (km/s) 

Fig. 1. Cross-correlation profile of VUMi obtained with 
the default KOIII (crosses) and M4V (diamonds) masks, 
during night N7 at phase ^0.99 (the night number refers 
to Table 2 of Paper III); the radial velocities as given 
by gaussian fits are: —171.2 km s _1 (K0, solid line) and 
-171.9 km s" 1 (M4, dashed line) 



Provence (France). The observations and the star sam- 
ple are described in detail in Paper III (Tables 1 and 2). 
The ELODIE spectrograph is designed to perform very ac- 
curate radial-velocity measurements by cross-correlating 
the stellar spectrum with numerical masks. In a first 
step, the cross-correlation functions (CCF) of the pro- 
gram stars were computed using two different numerical 
masks included in the ELODIE reduction software: (i) a 
KOIII mask constructed by Baranne et al. (1996); (ii) a 
M4V mask constructed by Delfosse et al. (1999) from an 
ELODIE spectrum of Barnard's star (Gl 699) applying 
the method of Baranne et al. (1979). 

2.1. Comparison of the CCFs obtained with the K0- 
and M4-templates 

The CCFs obtained for a given LPV star with the K0- and 
M4-templates are often quite different. Almost all situa- 
tions may be encountered. The most common are: (i) sin- 
gle peak with similar or different contrasts (Figs. [I] and g, 
respectively); (ii) double or asymmetrical peak with both 
templates (Figs. || and[|); (iii) double peak with the K0- 
template, red peak with the M4-template (Fig. ||); (iv) 
blue peak with the KO-template, double peak with the 
M4-template (Fig. ||); (v) doubtful single peak (Fig. fj]), 
doubtful double peak (Fig. |sj> or noisy profile (Fig. ^) with 
the KO-template, single peak with the M4-templatc. 

2.2. Interpretation: a first step towards tomography 

The variety of CCF pairs obt ained for a given star with the 
K0- and M4-templates (Sect. 2.1 ) may be easily explained 
if one admits that the KO-template probes deeper layers 
than the M4-template (the validity of this assumption will 
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Fig. 2. Same as Fig. (j]for RUHer (night N9, phase 0.02). 
Radial velocities are: —25.0 km s _1 (K0, solid line) and 
-25.2 km s _1 (M4, dashed line) 
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Fig. 3. Same as Fig. | for RT Cyg (night N2, phase 1.10). 
Radial velocities are: —130.7/— 111.0 km s" 1 (K0, solid 
line) and -134.9/-111.6 km s" 1 (M4, dashed line) 
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Fig. 4. Same as Fig. | for UZHya (night N7, phase 1.09). 
Radial velocities are: +28.6/+44.8 km s" 1 (K0, solid line) 
and +27.4/+41.0 km s _1 (M4, dashed line) 
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Fig. 5. Same as Fig. | for VTau (night N15, phase 2.01). 
Radial velocities are: +56.0/+77.8 km s _1 (K0, solid line) 
and +73.8 km s -1 (M4, dashed line) 
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Fig. 6. 

Radial 
-249.5 



Same as Fig. [j] for RYCep (night N4, phase 1.03:). 
velocities are: —245.1 km s _1 (K0, solid line) and 
/-229.3 km s" 1 (M4, dashed line) 
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Fig. 7. Same as Fig. for RPeg (night N15, phase 2.00). 
Radial velocity is: +21.5 km s _1 (M4, dashed line) 
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Fig. 8. Same as Fig. for TCep (night Nl, phase 0.23). 
Radial velocity is: —10.9 km s _1 (M4, dashed line) 
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Fig. 9. Same as Fig. gfor RTDra (night N12, phase 1.13). 
Radial velocity is: —104.0 km s _1 (M4, dashed line) 



in fact be demonstrated in Sect. 3.3.2 ) . As may be seen on 
Fig. [l(], the two templates catch different groups of spec- 
tral lines. Depending on the phase and on the spectral 
and chemical types, both, none or only one of the default 
templates will be able to detect the velocity discontinu- 
ity between the deeper, ascending layers and the upper, 
inf ailing layers. 

It will be shown in Paper III (Table 3 and Sect. 4.3.1) 
that double peaks are observed more often with the K0- 
template, whereas the M4-template often reveals only a 
single red peak (as in Fig. |J). This conclusion is a natural 
consequence of the fact that (i) the M4-template is more 
sensitive to the cooler, outermost layers which generally 
exhibit an infall motion around maximum light (hence the 
corresponding CCF is a red peak), (ii) the KO-template 
is more sensitive to hotter, i.e. inner, layers and, hence, 
is more likely to catch the discontinuity between the as- 
cending and infalling layers (corresponding to a double- 
peak CCF). A similar conclusion was reached by Hinklc 
et al. (1984), since the CO Av = 3 lines 1.6 /jm, forming 



4 



R. Alvarez et al.: Envelope tomography of LPV stars. II 



6.0x10 - 






5700 5705 
I ■: 1 1 I-]] 



Fig. 10. Part of the spectrum of /iGem (M3), observed 
during night N6. The thick segments represent the "holes" 
of the K0- and M4-templates 

deeper in the photospere than the blue- violet optical lines, 
exhibit line doubling while the latter do not. 

Sometimes however, the M4-template is also able to 
detect a double peak. In these cases, the KO-template may 
reveal either a double peak as well, the blue peak being 
then more pronounced than with the M4-template (Figs. || 
and|^), or solely the blue peak (Fig. ^|). 

When the star is at minimum light, the discontinuity 
front has left the photosphere so that the line doubling 
is no more visible. The late spectral types exhibited by 
most LPVs at minimum light prevent us from obtaining 
a well-defined CCF with the KO-template, which becomes 
useless. On the contrary, a single-peak CCF can still be 
obtained with the M4-template (Figs. || and ||) which thus 
delivers a well-defined CCF over the complete light cycle 
of (oxygen-rich) LPVs. 

Another interesting conclusion may be drawn from the 
comparison of the CCFs obtained with the K0- and M4- 
templates for stars which exhibit double peaks: the double 
peak appears at later phases with the M4-template than 
with the KO-template. As Fig. 
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indicates, a weak sec- 
ondary peak appears in RTCyg at phase 0.10 when ob- 
served with the M4-template as compared to phase 0.88 
with the KO-template. The same conclusion is reached for 
RYCep, a Ke-MOe Mira (Fig. fj]j). 

This phase lag can be interpreted once again in terms 
of the Schwarzschild scenario: since the shock wave prop- 
agates outwards and the M4-template scans higher layers 
than the KO-template, such a phase shift is not surpris- 
ing. It explains why the doubling is often observed after 
maximum light with the M4-template. 

For the same reason, the red peaks observed with the 
M4-template are almost systematically blue-shifted by 2- 
4 km s _1 as compared to the ones observed with the KO- 
template (see Figs. ^ and ||, and Table 3 of Paper III). This 
difference may be interpreted in terms of a velocity gradi- 
ent: since the infalling layers scanned by the M4-template 
are located above the ones scanned by the KO-template, it 
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Fig. 11. Sequence of CCFs for RT Cyg obtained with the 
KOIII mask (left side) and the M4 V mask (right side) dur- 
ing a two-month- long monitoring in August/September 
1999. The labels denote the visual phases 



seems plausible that these distinct layers fall with different 
velocities, reflecting the acceleration due to the gravity. 

The results discussed in this section thus show that the 
simultaneous use of the K0- and the M4-templates pro- 
vides insights into the velocity field of LPV atmospheres. 
Even more interestingly, the capability of the different 
templates to probe different layers in the atmosphere, as 
determined by their different spectral (i.e. temperature) 
adequacy, provides a hint on how to develop a method to 
investigate the propagation of the shock wave in LPV at- 
mospheres. This method, which consists in designing spe- 
cific numerical templates to perform the "tomography" of 
the atmosphere, is described in detail in Sect. |^. 

3. The tomographic technique 

The basic idea is to design a set of numerical templates 
that probe different atmospheric layers, in order to follow 
the propagation of the shock wave through the atmosphere 
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Fig. 12. Sequence of CCFs for RYCep with the default 
KOIII mask (left side) and the M4V mask (right side). 
The labels denote the visual phases. They are ordered ac- 
cording to increasing fractional phase within each cycle, 
and are centered around maximum light. The vertical lines 
are a guide to the eye 



(a similar technique was applied to BA-type supergiants 
by Kaufer et al. 1997). This technique takes advantage 
of the ability of the cross-correlation method to reveal 
complex profiles despite the severe crowding of the opti- 
cal spectra. The selective sensitivity of the default K0- 
and M4-templates demonstrates the potentiality of the 
method (Sect. ||). The crucial point of the tomographic 
technique is the selection of adequate sets of lines forming 
at a given depth in the atmosphere. 



3.1. Synthetic spectra 

The tomographic method rests on our ability to construct 
reliable synthetic spectra of late-type giant stars and to se- 




6000 6500 



Fig. 13. Synthetic spectrum with T cS =3500 K, log g=0.9, 
M=1.5 Mq. The dotted line is part of an optical spectrum 
of HD 123657, an M5-giant, from Serote-Roos et al. (1996) 



lcct the right set of spectral lines forming at a given depth. 
We compute synthetic spectra from static models of red 
giant stars in spherical symmetry (Plez et al. 1992; Plcz 
1992) in the spectral range 3850-6900 A. The synthetic 
spectra are computed at a high resolution (AA = 0.03 A, 
which corresponds to the binning step of the ELODIE 
spectra) using a set of routines especially designed for 
cool star atmospheric conditions ( Turbo spectrum] see Plez 
et al. 1993). The spectral synthesis includes the updated 
line data for TiO (Plez 1998), H 2 and VO as detailed 
in Bcssell et al. (1998) and Alvarez & Plez (1998), plus 
CH, C 2 , CN, NH, OH, MgH, SiH, CaH, ZrO, and atomic 
lines (Kurucz & Bell 1995). It results in more than 2. 10 7 
lines. The spectra are computed with CNO abundances 
and 12 C/ 13 C, 14 N/ 15 N and 16 0/ 18 isotopic ratios (re- 
spectively 13, 1250 and 500) typical of red giants (Smith 
& Lambert 1990). 

Figure [l3| shows an example of a synthetic spectrum 
(T c fF =3500 K). The dotted line is part of an optical spec- 
trum of HD 123657, an M5-giant, obtained by Serote-Roos 
et al. (1996) at a resolution of 1.25 A. The spectral region 
covered starts at 4800 A. 

3.2. Average depth of formation of spectral lines 



The static models described in Sect. 3.1 provide us with 



the monochromatic optical depths t\ as a function of ge- 
ometrical depth, the latter being actually represented in 
terms of x=logro, where To is an optical depth at some 
reference wavelength (here Ao = 1.2 /im). In static models, 
the monochromatic optical depth at any reference wave- 
length is indeed a monotonic function of the geometrical 
depth in the atmosphere. 

In the framework of the Eddington-Barbier approxi- 
mation, the emergent flux at wavelength A is supposed 
to come from the layer with t\ = 2/3. In the first im- 
plementation of the tomographic method presented here, 
we will assume that the line depression forms as well at 
that same depth. This assumption is valid for sufficiently 
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strong lines, as shown by Magain (1986). Using the expres- 
sion of the contribution function (CF)Q to the line depres- 
sion in intensity, Cr, Magain (1986; see also Gurtovenko 
et al. 1974) has however stressed that the region where the 
line depression is formed may be different from the region 
of origin of the emergent intensity, as described by the 
CF to the line intensity Cj (see Albrow & Cottrell 1996 
for the generalization to the line depression in flux, in a 
non-static atmosphere). Magain demonstrates that Cr is 
markedly different from Cj in the case of faint lines formed 
in classical stellar atmospheres (his Fig. 1). In the case of 
a stronger line Cr and Cj are more similar, so that the 
line flux and depression originate in the same layers (his 
Fig. 2). Cr and Cj are always close to each other in layers 
above r c = 10~ 2 (his Fig. 3), i.e. for sufficiently strong 
lines. Thus masks probing the deepest layers are the m ore 
likely to yield poor results, as indicated in Sect. 3.3.2 (see 
also Fig.gol). 



3851.5 
A (A) 



Fig. 14. Illustration of the design of a synthetic template. 
The present template (thick dashes at the bottom of the 
plot) scans the lines forming at depths x such that — 1.3 < 
x < —3.3 (see text) 



3.3. Designing synthetic numerical templates 

The synthetic templates are built in two steps. A complete 
spectrum, with all atomic and molecular lines included, is 
computed. The monochromatic optical depth t\ is cal- 
culated at each depth of the atmospheric model, for each 
sampled wavelength. The different masks Mj are then con- 
structed from the collection of N wavelengths Xij (1 < 
j < N) such that Xi < x(r\ i . = 2/3) < Xi+i — Xi + Ax, 
where x=\ogTo is some reference optical depth used to 
parametrize the geometrical depth. Each mask Mj thus 
probes lines forming at depths in the range X4, Xi + Ax. 

Figure [l4| illustrates how the synthetic templates are 
built: for the sake of clarity, only a very small portion 
(3 A) of the depth function x = x{t\ — 2/3) is represented 
and the constant Aa; has been assigned a very large value 
[Ax=2). The depth function provides the reference optical 
depth To = 1CP - or, equivalently, the geometrical depth 
- at which the monochromatic optical depth reaches 2/3. 
The depth function x — x(t\ — 2/3) therefore identifies 
the layer from which the monochromatic flux at wave- 
length Xij emerges. The triangles in Fig. ^ denote the 
local minima of the depth function (they correspond to 
the spectral lines): when these minima fall in the range 
scanned by a given template (i.e., between the two hori- 
zontal dashed lines), a corresponding "hole" is created in 
this template (thick dashes at the bottom of the plot). 

The constant Aa; was chosen after evaluating the con- 
tribution function to the emergent flux for some specific 
atomic lines. It appeared that the typical width Wcf of 
the CFs expressed in the a>scale is about Wcf=1-5. The 
Nyquist-Shannon theorem of elementary signal theory in- 
dicates that the optimal sampling of a given function is 
obtained for Ay/2 where Ay is the "resolution" of the 

1 The contribution function gives the relative contribution of 
the different atmospheric layers to an observed quantity which 
can be the emergent intensity or flux or the line depression in 
the surface intensity or flux 
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Fig. 15. Number of lines as a function of their formation 
depth (expressed in the log To scale; To = 1.2/zm) for a 
model at T eff =3500 K 

signal. A value of Ax=Wcf/2=0.75 was thus adopted to 
optimize the thickness of the atmospheric region scanned 
by each template Mj. 

Three sets of numerical templates were built from 
models with T cff =2800, 3500 and 4250 K. Figure |l| shows 
the distribution of the depths of line formation for a 
model with T o fT=3500 K: as expected, the number of 
lines increases with the atmospheric depth (as long as 
log to < —2). The lines forming in the outer atmospheric 
layers are by far less numerous, but they are also more 
intense. This justifies to keep the same Ax for all masks 
Mf. the number of lines and their intensities in any given 
mask vary in opposite directions so as to provide simi- 
lar signal-to-noise ratios for the CCFs obtained with the 
different masks (see Figs. [ll| [l?] and |l8| ). 

When these masks constructed with the complete set 
of lines were applied to non-variable stars, spurious sec- 
ondary peaks appeared in the CCFs obtained with some of 
the masks. These spurious peaks originate in inaccuracies 
in the wavelengths of some molecular transitions, that are 
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difficult to correct. Therefore, the decision was taken to 
rely solely on the supposedly more accurate atomic lines 
for building the templates. A new set of templates M[ 
were designed following the method described above, ex- 
cept that they were based on synthetic spectra includ- 
ing only atomic lines. Since important molecular opacity 
sources are lacking in these synthetic spectra, the asso- 
ciated depth functions may no longer be reliable, casting 
doubts on the line assignments leading to the new set of 
masks M[. 

To circumvent that difficulty, the final masks Ci 
are constructed by keeping all the holes from the M[ 
templates that appeared as well in the Mi templates. 
This way of doing guarantees that only the holes 
associated with atomic lines are kept (as they are 
present in the M[ masks), and that they are located 
in the right template, i.e. the template which scans 
the right depth range [xj,Xj+Aa;] (as they are also 
present in the Mj masks). This two-step procedure has 
the advantage to be easy to manage and avoids the 
enormous task of directly identifying which line (among 
the millions) contributes to a given hole before deciding 
to keep it or not. Table |l| gives for each atmospheric 
model the total range of reference optical depths x 
scanned by the masks, the depth step Ax and the num- 
ber of lines per mask Ci. The masks d are available at 



We defer to a forthcoming paper the test of how well this 
assumption is satisfied in practice. At this stage, we can 
only provide empirical, a posteriori evidence that this con- 
dition does not seem to be very badly violated at least. 
This empirical evidence relies on the fact that (i) the CCF 
shapes vary gently and smoothly from one mask to the 
next (see Figs. |l7] and |l^), and (ii) the evolution of the 
CCF shape from one mask to the next reflects what is 
expected in the framework of the Schwarzschild scenario 
(see Fig. 4 of Paper I). 



3.3.1. Application to non-variable stars 

In order to check whether our synthetic templates are reli- 
able, the CCFs of the non- variable stars of the sample have 
been computed with the newly designed masks. The se- 
quence of CCFs obtained for fi Gem, a M3 giant, with the 
3500 K templates is presented in Fig. 16, The CCFs ob- 
tained with the K0- and M4-templates are also displayed 
for comparison. 

It can be seen that the synthetic masks Ci yield the 
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At this point, a comment should be made about the 
validity of the use of static models to design templates to 
be used on dynamical atmospheres. Let us first remark 
that, if one aims at uncovering the complex line profiles 
of LPVs caused by the dynamics of their envelope, the 
correlation of the observed spectrum with a mask con- 
structed from a dynamical model is not the way to go. 
In the ideal case where the dynamical model used to con- 
struct the template would perfectly match the real star, 
the resulting CCF would be very sharp and single-peaked. 
To reveal the various velocities characterizing the different 
layers moving in the envelope, one needs to start instead 
from templates constructed from static models! 

The real difficulty in the application of 'static' tem- 
plates to dynamical atmospheres is elsewhere. It is clear 
that the contribution functions to the line depression may 
acquire several local maxima in a dynamical atmosphere 
(see e.g., Albrow & Cottrell 1996). This is precisely why 
the line shape acquires complex profiles in dynamical at- 
mospheres, which simply reflect the complex shape of the 
contribution function. The condition that needs to be sat- 
isfied for the tomographic method to yield reliable results 
with 'static' masks is the following: the contribution func- 
tions of all the lines that were assigned to a given static 
mask should change in the same way when passing from 
a static to a dynamic atmosphere. In other words, the 
lines supposed to form in the same layers in a static at- 
mosphere should still form in the same (possibly becoming 
non-connex) layers in a dynamical atmosphere. Or equiv- 
alently, there should be no swap of lines between masks 
when passing from a static to a dynamical atmosphere. 



same radial velocity to within 0.45 km s r.m.s. (exclud- 
ing the mask C8 which delivers a very broad CCF). Fig. [l^ 
also shows that the contrast of the CCFs increases along 
the sequence C\ to C%. This is not surprising since, from 
apffy.to Cs, each template scans stronger and stronger lines 



(i.e. forming closer and closer to the surface), and the CCF 
contrast reflects the average intensity of the group of lines 
scanned. It is also remarkable to note that the "signal- 
to-noise" ratio of the CCFs (ratio of the fluctuations of 
the continuum level to the depth of the peak) remains 
more or less constant. Also note how the wings of the 
CCF strengthen from C\ to Cs, revealing a similar - and 
expected - change in the profiles of the lines scanned. 

A further confirmation of the validity of our method 



is provided a posteriori by the discussion of Sect. 3.3.3 



which shows that the CCF correctly reproduces the com- 
plex shapes observed in a dynamical atmosphere for the 
lines probed by the considered mask (see Fig. |2l|). 



3.3.2. Application to variable stars 

Now that the reliability of the synthetic masks has 
been assessed from their application to non-variable stars 
(Sect. 



3.3.1 ), they may be applied to variable stars. 
Figures i"|l7| and |j~g| present the sequence of CCFs obtained 
with the synthetic templates of the T e ff=3500 K series 
for Z Oph and VTau during night Nl (phase 0.08) and 
N4 (phase 0.08) respectively; they clearly reveal a tran- 
sition from a single blue peak (templates C1-C3) to a 
red one (Cs) through double profiles (C4-C7). Again, the 
Schwarzschild scenario reveals itself, but this time in a spa- 
tial rather than temporal sequence. The deeper the mask 
scans, the more intense the blue peak is (and conversely 
for the red peak). With the synthetic templates, it is thus 
possible to follow the change of the CCFs from one layer 
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Table 1. The synthetic templates 
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Fig. 16. Sequence of CCFs obtained for /iGem with the synthetic templates of the T e g=3500 K series. Labels refer to 
the template names (see Table [[]). The sequence goes from the innermost (top left) to the outermost (bottom right) 
layers. The CCFs obtained with the K0- and the M4-templates are also displayed. Note how the CCF wings and 
contrast increase from C% to C$ 



to the next, i.e. to monitor the velocity field across the 
atmosphere. 

It is also interesting to note that the CCF obtained 
with the ^-template (i.e. deep layers) is very similar to 
the one obtained with the KO-template, whereas the Cq- 
and C7-templates (outermost layers) give profiles similar 
to those obtained with the M4-template. This similar- 
ity corroborates the conclusion of Sect. ^ on the different 
depths probed by the two default K0- and M4-templates. 

Fig. [li] provides another illustration of the efficiency 
of the synthetic templates: while each of the default tem- 



plates yields a single peak (although rather asymmetrical 
for the KO-template) for RT Cyg at phase 1.23, the C3- 
and ^-templates clearly reveal that line-doubling is actu- 
ally present in that star for the specific set of lines probed 
by the masks C3 and C4. We refer the reader to Paper I for 
the application of the tomographic masks to a sequence 
of spectra of RT Cyg covering phases 0.81 to 1.16, which 
clearly reveals the outward propagation of the shock front. 

The tomographic masks yield somewhat poorer results 
for late-type LPVs. As an illustration, Fig. ^ shows the 
application of the tomographic masks to the late-type 
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Fig. 17. Sequence of CCFs obtained for Z Oph (night Nl, phase 0.08) with the synthetic templates of the T e g=3500 K 
series. Labels have the same meaning as in Fig. |l^ 



LPV RU Her (M6-M9). Like the M4V mask, the outer 
masks C4 to C8 (of the 3500 K set) yield a well-contrasted, 
but sometimes quite asymmetric, CCF. The innermost 
masks CI and C2 provide a CCF that is, however, totally 
useless. The C3 mask presents several secondary peaks 
(that are reminiscent of those obtained with the KOIII 
template). The peak at V ~ 50 km s^ 1 remains visible in 
the C4 and C5 CCFs, but it is difficult to assess whether 
it traces a real layer or whether it must be ascribed to 
some kind of 'correlation noise' due to the mismatch be- 
tween the template and the observed spectrum. A future 
implementation of the tomographic method (lifting some 
of the simplifying hypotheses described in Sect. |^) should 
hopefully improve the results for late- type LP Vs. 



remarkably well the shape of the A 6358.69 line. On the 
contrary, the shapes of the CCFs obtained with the K0- 
and the M4-templates are noticeably different, since they 
were shown to correspond instead to the C2 and Cq — CV 
templates respectively (Sect. |3.3.2| ). This is clear evidence 
that the CCFs obtained with the synthetic templates cor- 
rectly reflect the shape of the set of lines they scan (at least 
in the case of early-type LPVs). This result may be seen 
as an a posteriori validation of the tomographic method 



as a whole (at least for early-type LPVs; see Sect. 3.3.2 
and Fig. |20| for a discussion of late-type LPVs) and of 
the various assumptions involved (use of 'static' masks, 
Eddington-Barbier approximation for the depth of forma- 
tion of spectral lines...) as discussed in Sect. [| 



3.3.3. Individual lines and cross-correlation profiles 

A definite proof that the complex CCF profiles of LPVs as 
observed for instance on Figs. or H are not an artefact 
of the method, is provided by the comparison of the CCF 
profile with the shapes of individual lines probed by the 
corresponding mask. For this purpose, an early-type LPV 
must be selected, to avoid too severe spectral crowding. 
Z Oph (K3 at maximum) has the earliest spectral type 
in the sample. Fig. |l] compares the A 6358.69 Fe line at 
4 different phases (bottom row) with the CCFs obtained 
with the K0-, M4- and Cs-templates (from top to bottom), 
since the A 6358.69 line belongs to the Cs-template. The 
CCFs obtained with this template (third row) reproduce 



3.3.4. Comparison of CCFs obtained with synthetic 
templates constructed from models of different 
temperatures 

This section briefly addresses the question of the influ- 
ence of the model used to construct the synthetic tem- 
plates on the resulting CCFs of LPVs. Browsing the large 
database containing the CCFs computed with the tomo- 
graphic masks of the T c g = 4250, 3500 and 2800 K series 
for all 315 available LPV spectra (Table 3 of Paper III) 
has shown that (i) the masks based on the 4250 K and 
3500 K series are the most appropriate to perform the to- 
mography of LPVs (since the CCFs obtained with masks 
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Fig. 18. Sequence of CCFs obtained for VTau (night N4, phase 0.08) with the synthetic templates of the T e gf=3500 K 
series. Labels have the same meaning as in Fig. |l^ 



Ci, C2 of the 28 00 K series are generally noisyQ- see also 
Fig. ^ in Sect. 3.3.2| ); (ii) no information that is not al- 
ready present in the CCFs of the 4250 K series is obtained 
with the CCFs from the other two series. Schematically, 
the CCFs of the three series match each other after trans- 
lating one sequence with respect to the other, i.e., the CCF 
of template (7^(4250 K) is almost identical to the CCF of 
template Ci_i(3500 K). 

4. Summary 

A tomographic method has been described that makes 
it possible to follow the propagation of the shock wave 
across the photosphere of LPV stars. The method relies on 
the correlation of the observed spectrum with numerical 
masks probing layers of different atmospheric depths. The 
tomographic masks are constructed from synthetic spec- 
tra of red giant stars, that provide the depth of formation 
of spectral lines. When applied to Mira stars around maxi- 
mum light, they reveal that the deeper layers are generally 
characterized by blueshifted absorption lines (translating 
their upward motion) , whereas the uppermost layers gen- 
erally exhibit redshifted absorption lines (translating their 
infalling motion). A double-peak profile is then found in 
intermediate layers, where the shock front is located. At 



note that this poor result may be due to the condition 
t\ = 2/3 used to define the depth of formation of the line 
depression, which was shown in Sect. 3.2 to be poorly satisfied 
in the case of faint lines 



later phases, the shock front is seen moving towards upper 
layers, until it leaves the photosphere. 

A number of checks validating the tomographic 
method have been presented: (i) the tomographic masks 
yield single-peaked CCFs at the same radial velocity (with 
a r.m.s. deviation of about 0.5 km s _1 ) when applied to 
non-pulsating stars; (ii) the default K0- and M4-templates 
behave approximately like the C2 and Cq — C7 templates, 
respectively, of the T e g — 3500 K series; (iii) when checked 
on a spectrum of an early-type LPV that is not too 
crowded, the CCF profile correctly reproduces the com- 
plex shape of the lines belonging to the corresponding 
mask. 

The application of the tomographic method to a se- 
quence of spectra around maximum light makes it possi- 
ble to follow the outward propagation of the shock front 
(this was demonstrated in Paper I for RT Cyg). 

The method offers interesting perspectives for the 
study of the dynamics of LPV atmospheres. However, 
in its current state, its quantitative predictions (concern- 
ing for example the velocity of the shock front) remain 
rather limited, mainly because the geometrical radius as- 
sociated with each mask is currently unknown. Spectro- 
intcrfcrometric observations would be necessary to make 
progress in that direction. Such measurements will prob- 
ably become feasible with the new generation of interfer- 
ometers that will be available in the near future. 
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Fig. 20. Sequence of CCFs obtained for the late-type LPV RUHer (M6-M9; night N10, phase 0.08) with the synthetic 
templates of the T e ff=3500 K series. Labels have the same meaning as in Fig. [ji] 
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Fig. 21. Comparison of the A 6358.69 Fe line as seen in the spectrum of Z Oph at phases 0.08, 0.17, 1.04 and 1.21 
(bottom row) with the CCFs obtained with the K0-, M4- and CVtemplates (from top to bottom) 



